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Graphene p-n junctions with nonuniform Rashba spin-orbit coupling 
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Linear conductance of graphenc-based p-n junctions with Rashba spin-orbit coupling is considered theoret- 
ically. A square potential step is used to model the junctions, while the coupling is introduced in terms of 
the Kane-Mele model (C.L. Kane and E.J. Mele, Phys.Rev.Lett. 95, 226801(2005)). The main objective is 
a description of electronic transport in junctions where Rashba parameter is nonuniform. Such a nonunifor- 
mity can appear when graphene is asymmetrically covered with atomic layers, or when Rashba coupling is 
strongly dependent on electric field. It is shown that conductance is significantly modified by the considered 
nonuniformity, which is most clearly manifested by an anomalous minimum at a certain potential step height. 



Owing to the electric field effect 1 , various all-graphcnc- 
based heterostructures can be easily made by means 
of electrostatic gates. The simplest structures of this 
kind are p-n junctions, which have been already anal- 
ysed theoretically 2 and investigated experimentally 3 . 
Generally, transport characteristics of such junctions 
significantly depend on the spin-orbit interactions in 
graphene 4-6 . This is mainly due to the fact that ac- 
cording to the commonly used Kane-Mele model 7 , the 
intrinsic coupling opens a band gap, while the Rashba 
coupling (RSOC) splits both the conduction and valence 
bands. In this paper we focus on the latter, motivated by 
recent results showing that it can be made significantly 
dominant in graphene, as pointed out in Refs. 12-18, 
where RSOC parameters of the order of 10 meV were 
reported. 

The parameter of RSOC has been so far treated in 
transport investigations as a constant throughout the 
structure and also independent of the gate voltages used 
to make p-n junctions 4 ' 5 . However, this parameter can be 
generally nonuniform, which may occur when graphene 
is covered with an atomic layer from bottom in one part 
of the junction and on top in the second part 12-18 , or 
when RSOC is strongly dependent on electric field 8-10 . 
This problem is addressed in the present paper, where 
we show that conductance of such p-n junctions differs 
from that of junctions with constant spin-orbit coupling. 
We focus especially on a situation in which the RSOC 
parameter changes sign at the border between the two 
parts of the junction. All this gives rise to anomalous 
behaviour of the corresponding conductance. 

In order to describe charge carriers in graphene we 
employ the effective low-energy Hamiltonian , Hq = 
—iTivp(a x d x + o-ydy) <g> So, where vp is the electron ve- 
locity. We use the notation with cr a and s a being the 
Pauli matrices (for a — x,y,z) in the pscudospin and 
spin spaces, respectively, and <jq and sq denoting the 
corresponding unit matrices. We restrict our consider- 
ations to a single electronic valley, neglecting thus any 
inter-valley scattering. RSOC is taken in the form intro- 
duced by Kane and Mele 7 , Hr = \(a x ® s y — a y ® s x ), 
where A is the coupling parameter. The full Hamiltonian, 



H = Hq + Hr, has the following eigenvalues: 



E = IX + s^/X 2 + (fc 2 + k 2 y )v 2 F h 2 



(1) 



where k is the two-dimensional wave vector, s = 1 (s = 
— 1) for the conduction (valence) band, and I = ±1 is 
used to distinguish between the two subbands. 

The p-n junction is described by the term, V(x, y) = 
Vo<d(x)<r <E> Sq, to be added to the full Hamiltonian. Let 
us denote the Fermi energy (measured from the neutrality 
point in the left part of the junction) by Ep. 

Assume now a particle of energy E at the Fermi level, 
E = Ep, in the subband I is incident on the poten- 
tial step at an arbitrary angle 4>. The corresponding 
wave function has the form of a plane wave multiplied 
by a column vector derived in Rcf. 4, ip E ^ t (x, y) = 
\sk x ,E,4>,l) e^p[isk x x\exp[ik y y] = ip^^x) exp[ik y y], 
where k x is a positive solution of Eq. (1) (note s ensures 
that group velocity points in the right direction). In the 
same way one can build the wave functions for carriers re- 
flected and transmitted into the subband I; riip E , ^x, y) 
and tiip E _ Vo ^ i(x, y), and for carriers reflected and trans- 
mitted into the other subband, I' ^ Z; r[%\i E ^ v (x,y) and 

t'l^E-Vo 4> l'( x > v) m Here, n, U, r[ and t\ are the corre- 
sponding reflection and transmission amplitudes. Owing 
to the translational symmetry along the y direction, the 
k y component of the wave vector is conserved. 

To find the linear conductance we need to know the 
amplitudes ti and t[ , which can be found from the conti- 
nuity condition at the boundary x = 0, 

Vl W (z = 0) + M^^x = 0) + r[xl) EM ,{x = 0) 



W-vw,i(a: = o) + tlv 



E-V ,ef>,l' 



(x = 0). 



Having found the transmission amplitudes and taking 
also into account the respective group velocities, one ob- 
tains the corresponding transmission probabilities TJ and 
T[. The full linear conductance can be then obtained by 
integrating over all incidence angles, 



G 



2k F e 2 W 



x ^ k 
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F 



r/)cos(«. (3) 
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FIG. 1. Band structure of graphene for |Ai| = 10 meV, a -1 = 
0, Ef — 50 meV, and Vo = 100 meV. Left part shows the 
spectrum on the incident side of the junction (where Ai = 
10 meV), while the right part shows the spectrum on the 
transmitted side (where Ai = —10 meV). The solid (dotted) 
lines distinguish subband with opposite index /. 



In the following, the conductance will be normalized to 
the factor Go defined as Go = 2e 2 Wkp/(hir), where Uf 

is an average of incident Fermi wave vectors kp , while e 
is the electron charge and W is the sample width. 

The RSOC parameter A can in general be separated 
into two terms; A = Ai + \%{E). The first term is inde- 
pendent of electric field - it may be due to the adjacent 
atomic layers (substrate/cover layers). The second term 
can be controlled by external electric field E (gate volt- 
ages), X2(E) = ol\Z ■ E, where z is a unit vector perpen- 
dicular to the graphene plane, and ai is a relevant pa- 
rameter. The Fermi level in graphene is also controlled 
by electric field, and this dependence can be expressed 
as Ef = sign(z • E)y/\z ■ E\a.2, where a.2 is a parameter, 
while Ef = Ef for x < and Ef = Ef — Vo for x > 0. 
Thus, when the Fermi level is tuned by the gate voltage, 
the RSOC parameter A has to be adjusted according to 
the formula 



A = Ai + sign(E F ) 



-1 fp2 



(4) 



where a = a^/a\. In the following we will consider two 
options for the sign reversal of the RSOC parameter at 
the potential step. First situation is when oT x = and 
Ai is positive on one side of the step and negative on 
the other side. The second possibility appears when a -1 
is sufficiently large, while Ai is small (we take Ai = 
for clarity of the argument) and constant through the 
structure. We begin with the former. 

In Fig. 1 we show the band structure of graphene on 
the two sides of the junction for |Ai| = 10 meV. The 
assumed value is typical for graphene covered by atomic 
planes 12-18 . The parameter Ai is positive on the incident 
side and negative on the transmitted one (an opposite 
case would be fully symmetrical). Note, that apart from 
an up-shift, the spectrum on the transmitted side is in- 
verted with respect to that on the incident side. This 
inversion is due to the sign reversal of the RSOC param- 
eter Ai at the potential step. 

Figure 2 shows how the corresponding conductance G 
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FIG. 2. Normalized conductance versus potential step height 
Vo for | Ai| as indicated, with Ai positive on the incidence side 
and negative on the transmitted side. The other parameters 
are as in Fig.l 



depends on the potential step height Vq for four differ- 
ent values of |Ai|. The curve corresponding to Ai = 
serves as a reference. For the two curves correspond- 
ing to the lower values of |Ai| ^ the incident carri- 
ers come from both conduction subbands (as in the case 
shown in Fig. 1). Therefore the corresponding conduc- 
tance reaches values above G/Go = 1. There is a general 
trend in these curves that, when Vo moves from the global 
minimum at 50 meV towards larger values, the conduc- 
tance grows until it saturates, similarly to the |Ai| = 
case. When Vq moves towards lower values in turn, the 
conductance also grows at first and finally saturates, but 
the maximum present for Ai = OatVb = Ois sup- 
pressed. This reduction is due to band mismatch created 
by inverted band splitting on both sides of the junction. 
Moreover, in general Fermi level on the transmitted side 
crosses either one or two subbands. The point, where 
a new subband enters (or leaves) the conduction regime 
is associated with a kink in the conductance curve. The 
curve for |Ai| = 30 meV differs from the others. First of 
all it reaches values only up to G/Go = 1, since the car- 
riers on the incident side come from only one subband. 
Moreover, at Vq = the aforementioned band mismatch 
leads to a minimum equal zero. Furthermore, only one 
of the kinks mentioned above is visible in this case (at 
Vo = —10 meV). Finally, at Vo = 100 meV an anomalous 
conductance minimum occurs. 

Consider now the second situation, in which the RSOC 
parameter is dependent on electric field, and let us as- 
sume the simplest case where Ai = 0. In Fig. 3 we show 
the band structure of graphene on the incident and trans- 
mitted sides of the junction for indicated values of the 
step height Vo- Since the RSOC depends on the gate volt- 
ages (and therefore on the Fermi level and Vq), for each 
situation the RSOC parameter is adjusted according to 
the formula (4). This dependence can also lead to sign 
reversal of the RSOC parameter at the potential step. 
To emphasize this effect, we assumed a -1 significantly 
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FIG. 3. Band structure of graphene for Ai = 0, a = 0.25 eV, 
Ef = 50 meV, and for indicated values of the potential step 
height Vq. Part (a) shows the spectrum on the incident (left) 
side of the junction, while the parts (b) to (f) show spectra 
on the transmitted (right) side for Vb as indicated. 



larger than currently available data for graphene 10 . But 
one may expect, that this parameter can be increased by 
means of an appropriate substrate and/or cover layers. 

In Fig. 4 we present how the conductance G depends 
on the potential step height Vb for different values of 
a. The first three curves (corresponding to larger values 
of a) show similar behaviour, since in all these cases the 
incident carriers come from both conduction subbands. 
First, there is a minimum at Vq = 50 meV, correspond- 
ing to vanishing density of states on the transmitted side 
(RSOC vanishes as well, see Fig. 3(d)). When V in- 
creases from this value, the conductance grows until it 
reaches a maximum value larger than G/Go = 1, and 
then it suddenly drops. This drop occurs when one of 
the subbands on the transmitted side ceases taking part 
in conduction (transition from (e) to (f) in Fig. 3). When 
Vb changes from Vb = 50 meV towards lower values, the 
normalized conductance grows until it reaches a maxi- 
mum value, G/Go = 2, which is located at Vo = (no 
potential step, see (a) and (b) in Fig. 3). When Vb goes 
further to negative values, there is a drop in conductance 
of the same origin as the one discussed above. 

The curve for a = 0.08 eV in Fig. 4 corresponds to the 
case where the incident carriers come from one subband 
only. Therefore, the conductance is always equal to or 
smaller than Go, G/Gq < 1. As before, there is a max- 
imum at Vb = and a minimum at Vq = 50 meV. Ad- 
ditionally, an anomalous minimum appears at Vo = 100 
meV. This anomalous minimum is of the same nature as 
the one discussed in the case where a -1 = 0. 

In summary, we have considered the charge trans- 
port through p-n junctions in graphene with nonuniform 
RSOC parameter. The nonuniformity was either due to 
asymmetric location of adjacent atomic planes (on top in 
one part of the junction and at bottom in the other part), 




Step height V [meV] 

FIG. 4. Normalized conductance versus potential step height 
Vo for Ai = 0, Ef = 50 meV and indicated values of the 
parameter a. 

or due to a strong dependence of the RSOC parameter 
on electric field (and thus on the gate voltage used to 
generate carriers) . We have especially focused on the sit- 
uation when the RSOC parameter changes sign at the 
potential step. We have found a suppression of the con- 
ductance at characteristic points, which is a consequence 
of the inverted band splitting by RSOC. 
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